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Abstract 

The yields of dibaryon (r2r2)o+ in relativistic heavy ion collisions at three energies 
(AGS, SPS and RHIC energies) are calculated in the framework of statistical model. 
It's yield at full RHIC energy, y/s^N = 200 GeV, is predicted by means of extrap- 
olation. In the RHIC energy region (RHIC to full RHIC energies), it is in order of 
magnitude 10~^ — 10~^. The yields of hyperon and the ratios of (i7r2)o+ to Q.~ 
are also given. In the RHIC energy region, the yield of 0~ is in order of 10° - 10^ 
and the ratio of (r2r2)o+ to $7^ is ~ 1.0 x 10^^. 
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1 Introduction 

According to quantum chromodynamics (QCD) theory, there is the possibil- 
ity of the existence of some exotic multi-quark states, such as H dibaryon [1], 
and gluonic states. The growing interest of searching for such states lies on 
not only to examine the quantum chromodynamics (QCD) theory, but also 
to reveal the quark-gluon behavior in short distance and other new physical 
characteristics. But so far no convinced evidence is found in experiments to 
show their existence. Recently a new candidate, the strange dibaryon (nfi)o+, 
is suggested and studied based on the chiral SU(3) quark model [2-6]. This 
six-quark states with large strangeness is found to have more attractive char- 
acteristics than other six-quark states. First, it is a deeply bound state. Al- 
though the color magnetic interaction in the one gluon exchange term for this 
diomega system exhibits the repulsive feature, the large attraction stemming 
from the chiral gluon coupling and from the symmetry property of the sys- 
tem lead to a considerable binding. By the Resonating Group Method, the 
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binding energy of (r2f2)o+ is found to be as large as ^ 116 MeV, and the 
root-mean-square distance between two Q~s is 0.84 fm. Second, (QQ)o+ has 
a quite long lifetime of ~ 10~^° sec as it can undergo only weak decay, such 
as (f2r2)o+ 7r~ + + Q~ and (nO)o+ — >■ 7r° -|- S~ -|- Q~. As they are three- 
body decay, the final state phase-space would be suppressed. In the sudden 
approximation, the lifetime of (r]fi)o+ was found to be four times longer than 
the free lifetime of 0.822x10^^*^ sec. Apart from these conventional decay 
modes, the non-meson decay (r2f2)o+ — > S° -|- Q~ is also possible, the lifetime 
of (r2r2)o+ for this process is twice the free hfetime. Third, (nn)o+ has two 
negative charges. These properties make it easily be identified in experiment 
when it were created. 

Because of the large strangeness and very heavy mass, (f2r2)o+ is not likely 
to be produced in proton-proton collisions. While in relativistic heavy ion 
collisions, it is found that the strangeness production is enhanced and has been 
suggested as one of the possible signals of quark-gluon- plasma (QGP). Recent 
experiments by the WA97 Collaboration [7] and the NA49 Collaboration [8] 
showed the substantial enhancement of the (anti-)hyperon (A, S and Q) yields 
in 158 A GeV/c Pb-|-Pb central collisions relative to p+Pb collisions, and 
exhibited the enhancement pattern increasing with the strangeness content of 
the (anti-)hyperon. Therefore, the extreme environment of very hot and very 
dense block of matter provided by the relativistic heavy ion collisions may 
lead to the formation of exotic deeply bound objects composed of quarks or 
baryons with large strangeness, among them the dibaryon (f2r2)o+ is the most 
possible candidate. 

Theoretically, the yields of (QQ)o+ can be estimated by two ways: The kinetic 
simulation or microscopic transportation model, and the quantum statistical 
model or thermal model. Because the formation probability of (r2r2)o+ is very 
low and the cross sections related to (^2^2)0+ formations and decays are not 
well known, the calculations by kinetic simulation can not be well carried out. 
Contrarily, the statistical model can do the calculations easier. Based on this 
model, ref. 5 gave out an estimation of (Jin)o+ yields in Pb+Pb collisions 
at 158A GeV/c (SPS), which is in order of magnitude 10^^. In this paper, 
we will present its yields in other two reactions: Au-|-Au collision at 11. 6A 
GeV/c (AGS) and Au+Au collision at = 130 GeV (RHIC). And, based 

on the results of these three reactions, we try to predict the yield of (QQ)o+ 
for Au-|-Au collision at full RHIC energy, y^J/vw — 200 GeV, by means of 
extrapolation. In order to reduce the effect of system size to the extrapolated 
results, the above collisions are chosen to have the same or similar reaction 
system. The yields of Q~ and the ratios of (r2r2)o+ to Q~ are also given. 
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2 Thermal Production 



In relativistic heavy ion collisions, a large amount of particles including baryons, 
mesons, their resonance and anti-particles are produced and form a fireball of 
hadronic gases with high particle density and high energy density. The fireball 
undergoes the expansion until freeze out. The large amount of produced par- 
ticles and the violent collisions between them lead the gaseous particles inside 
the fireball to reach the chemical and thermal equilibrium. Such gas system 
can be described by the statistical model or thermal model. This model has 
been successfully applied in analysis of hadron yields and ratios at AGS and 
SPS energies [9-18], and recently the RHIC energy [19,20]. Here we use this 
model to estimate the yields of (r2r2)o+ in above three energies and predict it's 
yield in full RHIC energy, ^/snn — 200 GeV, by means of extrapolation. 

In the frame of grand canonical ensemble, the density of hadron species i reads 
as (in unit of c = h = ks = I) 



where T is the temperature of the fireball, s = yq"^ + rri^ is the energy of the 
particle, q and rrii are the momentum and mass, Qi is the spin-isospin degen- 
eracy. ± stands for the characteristics of statistics, -|- is for the Fermions and - 
is for the Bosons. The chemical potential is determined by //j = -Bj/ze + Si/is, 
Hb and are the chemical potentials of baryon and strangeness respectively, 
Bi and Si are the baryon charges and strangeness charges of the particle re- 
spectively. 

Since the temperature of the fireball is very high, it is reasonable to treat 
the hadron gas as the ideal gas obeying the Boltzmann distribution. In the 
condition oim/T ^ 1, the density of species i can be well approximated by 



We carried out the calculations for the following three relativistic heavy ion 
collisions: Au+Au collision at 11. 6A GeV/c (AGS), Pb+Pb colhsion at 158A 
GeV/c (SPS) and Au+Au colhsion at ^/i]^ = 130 GeV (RHIC). The ex- 
perimental data of the hadron yields and ratios are adopted from [10,12,19] 
and the references therein, respectively All hadrons including baryons, mesons 
and possible anti-particles with mass less than 2 GeV/c^ are included in the 
hadron gas of the fireball. The calculations are carried out within the frame 
of two-source statistical model (TSM) [17,18,20]. This model is developed on 




(1) 




P^= 9^(^f"e-^-^-^^)/^ . 



(2) 
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the basis of the standard statistical model, the single-source model (SM), 
and takes the non-homogeneous effect of the source into account. This effect 
is showed up in the experimental observation of low baryon components in 
the midrapidity region of relativistic heavy-ion collisions at energy of 158A 
GeV/c (see, e.g. [21]) and in the microscopic model calculations which show 
the central zone is more baryon dilute than outside zones [22,23]. This effect of 
"non-homogenous" implies the requirement of a multiple source model to bet- 
ter describe the source of a hadronic gas. The two-source model (TSM) is the 
simplest multiple source model. In TSM the whole source is roughly divided 
into two regions (sources): the inner source S2 and the outer source 5"!. The two 
sources arc supposed to reach equilibrium respectively and possess different 
temperatures, densities, chemical potentials and other thermodynamic char- 
acteristics. Especially the net strangeness density is no longer zero everywhere 
as in the single-source model but the total number of strangeness still keeps 
conserved. The two-source model significantly improves the agreement to the 
experimental data of hadron yields and ratios and exhibits a reasonable source 
structure [17,18,20]: At SPS and AGS energies, the whole source is composed 
of a small, hotter core (inner source) surrounded by a large, cooler halo (outer 
source). Most baryons are distributed outside, seems to be the projectile- and 
target-like components, while almost all anti-baryons are concentrated inside. 
Besides, in the two-source model approach, the average particle densities and 
the average energy density are effectively reduced and the strangeness is ef- 
fectively suppressed but without introducing extra corrections, such as the 
corrections of hard core excluded volume of the particles and the strangeness 
suppression. These corrections were generally used when the whole source is 
assumed to be a unique homogeneous source. 

By means of least square method and by fitting the experimental data, the 
main thermodynamic characteristics in the two sources, S2/S1, are obtained 
as follows: 

AGS energy the temperatures, T2/T1, are 144/91 MeV. The baryon chemical 
potentials, 1^32/ jJ'Bi, are 288/599 MeV. The strangeness chemical potentials, 
1^32/ ^si, are -70.0/116 MeV. 

SPS energy T2/T1=168/110 MeV, //b2/a«bi=14.4/387 MeV, /X52/a*5i=-14.8/52.2 
MeV. 

RHIC energy T2/T1=176/175 MeV, /iB2//^i3i=36.9/38.9 MeV, /i52//^si=10.0/10.7 
MeV. 

The Q~ yields in above three reactions are directly obtained through for- 
mula 2. Because of the large strangeness, -6, and very heavy mass, 3230 MeV, 
the (f2f2)o+ yields are very small. Therefore, it can be supposed that the ther- 
modynamic quantities would be little effected by the existence of (f2f2)o+. In 
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Fig. 1. Energy dependence of (J7J7)o+, ^~ and (J7J7)o+/n~ in relativistic heavy ion 
collisions. Points at fTill RHIC energy, y^ijv7v=200 GeV, are obtained by extrapola- 
tion. Points at ,/snn=75 GeV are obtained by interpolation. 

statistical model of grand canonical frame, all hadrons in the fireball includ- 
ing (r2^r2~)o+ and are in chemical and thermal equilibrium. According to 
formula 2, the yield ratio of (r2i7)o+ to has the relation of 

^ 9{nQ) . m^m \3/2 g-[{m(on)-M(nQ))-(mn-MQ)] ^3) 

9n run 

With formulas 2 and 3 the (r2f2)o+ yields can be obtained from the ^2" yields. 
The yields of (f2r2)o+ and Q~ together with the ratios are presented in Table 1 
(rows 2-4) and shown in Fig. 1. 



(put Table 1 here) 



Apparently, the yields of (Qil)o+ and Q increase steadily as the energy in- 
creases from AGS to SPS, and to RHIC. At RHIC energy, they are 5.8x10"^ 
and 5.8 respectively. The ratio, however, shows a saturation feature in RHIC 
energy region with value of ~ 1.0 x 10~^. 

With results in above three reactions, we can try to predict the results for 
Au+Au collisions at full RHIC energy, y^J/vw = 200 GeV, by means of ex- 
trapolation (with Lagrange formula). In order to reduce the effect of system 
size to the extrapolated results as possible, these three collisions have been 
chosen to possess the same or similar reaction system. The results are shown in 
Fig. 1 (last points) and also shown in Table 1 (last column). At RHIC and full 
RHIC energies, the yields of (r20)o+ are 5.8x10"^ and 1.2x10"^ respectively, 
well within the limits of the detectors used in RHIC energies. With the ex- 
cellent characteristics, (r2f2)o+, the marvelous particle, can be easily identified 
and detected when it is created. 
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3 Discussion and Summary 

In order to justify the reliability of the model predictions, we give out some 
results related to the strange particles and compare them to the experimental 
data or other model predictions. Table 2 and Fig. 2 show the yields of Q^, 
and their ratio At low energy the il^ yield is much lower than 

the yield, while fl increases faster than Q~ as the energy increases. So 
at somewhere they will have equal yield (in the figure, it is at ^/sNN ~ 150 
GeV). Table 3 and Fig. 3 show the ratios of the anti-strange particle to the 
strange particle. In Fig. 3 the results at AGS energy are for 11. 6 A GeV/c 
Au+Au collision used for extrapolations. While in Table 3 the results at AGS 
energy are for 14. 6A GeV/c Si+Au collision used to compare with the available 
experimental data[9]. For the strange baryons, the ratios possess the similar 
behavior as the energy increases but the magnitudes are different. Ratio for 
particle with more strange components is larger than that with less strange 
components. The ratio of kaon exhibits a different feature which saturates 
after = 70 GeV with value of 0.8. From the comparisons in Tables 2 

and 3, one can see that the model predictions are in very good agreement with 
the experimental data. 



(put Table 2 here) 



(put Table 3 here) 



Several points need be pointed out. Firstly, the yields and ratios obtained 
above are thermally produced, i.e. produced in a hot and dense fireball which 
is assumed to reach the equilibrium or local equilibrium. If taking the kinetic 
effect into account, some amount of (nn)o+ would be produced from two Q~s' 
collision, the total yield of (QQ)o+ will be a little increased. Secondly, the 
yields of Q~ predicted by statistical models generally lower than that given 
by experiments. So the yields of Q~ and (r2i7)o+ presented here are the lower 
limit. Thirdly, more importantly, the experimental data used in calculations 
for RHIC energy are taken from pseudorapidity interval I77I < 1 which is a 
small portion of the whole source of which the rapidity range is —5.4 < y < 
5.4. (This small portion of source, the "measured" source, is located in the 
central part of the source. In this region, the particles arc uniformly distributed 
and construct a single homogeneous source. So the two parts, SI and S2, 
in TSM are identical (see above the thermodynamic characteristics at RHIC 
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(put two figures in the same row) 



Fig. 2. Energy dependence of r2 , r2~ and f2 /O^ in relativistic heavy ion collisions. 
Points at full RHIC energy are obtained by extrapolation. Points at ■s/snn=75 GeV 
are obtained by interpolation. 

Fig. 3. Energy dependence of A/A, /Q^ and /K^ in relativistic 

heavy ion collisions. Points at full RHIC energy are obtained by extrapolation. 
Points at ■s/snn=75 GeV are obtained by interpolation. 

energy) and also identical to the source in the single-source model.) Therefore, 
if the data are taken from the whole source or a larger rapidity interval, the 
yields of Q~ will be raised and correspondingly, the yields of (QQ)o+ increases. 
For example, in a conservative estimation, suppose the yield of increases 
two times to reach 12, the yield of (f2f2)o+ increases the same times to reach 
1.2x10^'^. At full RHIC energy, their yields increase ~ 2.2 times and reach 
24 and 2.6x10^^ respectively. Obviously, a larger central rapidity range is 
favorable to measuring the interesting but rarely produced strange dibaryon, 
(Q0)o+. Of cause, the more exact results are expected after the data of full 
RHIC energy become available. 

So far experimental data are available only for three collisions at AGS, SPS 
and RHIC energies with the same or similar collision system. There is no data 
available within the large energy region between SPS and RHIC energies. The 
particle yields, not as the ratios, are quite sensitive to the "measured" source 
size or "measured" rapidity interval. Under these limitations, the yields are 
only roughly estimated and may have large deviation, but it is believable that 
the order of magnitude would not deviate far away. 

In summary, the yields of dibaryon (f2f2)o+, hyperon Q~ and their ratios in 
relativistic heavy ion collisions at three energies are calculated within the 
framework of two-source statistical model. The yields of (fin)o+ and Q' in- 
crease steadily as the energy increases, while their ratios display a saturation 
feature. By means of extrapolation, the yields and the ratio at full RHIC en- 
ergy are estimated. In the RHIC energy region, the (^2^2)0+ yield is in order 
of 10"^ - 10-3, the n- yield is in order of 10° - 10^ and the ratio of {nn)o+ 
to Q~ is ~ 1.0 X 10"^. And, to measure the rarely produced particles as the 
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strange dibaryon (QQ)o+, a larger central rapidity range is favorable. 
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FIGURE CAPTIONS 



Fig 1 Energy dependence of (r2r2)o+, fi^ and {QQ)q+ /fl^ in relativistic heavy 
ion collisions. Points at full RHIC energy, ^/snn— 200 GeV, are obtained by 
extrapolation. Points at GeV are obtained by interpolation. 

Fig 2 Energy dependence of fi"*", and fi^/i^" in relativistic heavy ion 
collisions. Points at full RHIC energy are obtained by extrapolation. Points 
at -s/snn—IB GeV are obtained by interpolation. 

Fig 3 Energy dependence of A/A, S^/S^, and K~ / in relativistic 

heavy ion collisions. Points at full RHIC energy are obtained by extrapolation. 
Points at y/sNN=75 GeV are obtained by interpolation. 
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Table 1 

Yields of (f2-r2-)o+ and n~, and their ratios in relativistic heavy ion collisions. 
Results at full RHIC energy, y^sjvjv=200 GeV, are obtained by extrapolation. 





AGS 


SPS 


RHIC 


full RHIC 


{n-n-)o+ 


9.18x10"^ 

0.0423 
2.17x10-5 


2.78x10-5 

0.444 
6.26x10-5 


5.78x10-'' 

5.81 
0.996x10-'' 


1.20x10-3 

10.7 
1.03x10-^ 



Table 2 

Yields of 0,~ and , and Jl^/fJ- in relativistic heavy ion collisions. Results at 
full RHIC energy are obtained by extrapolation. 





AGS 


SPS 


RHIC 


full RHIC 




4.23x10-2 


0.444 


5.81 


10.7 


n+ 


3.74x10-3 


0.166 


5.37 


12.0 




8.84x10-2 


0.378 


0.926 


1.04 






exp: 0.383±0.081[24,12] 


(0.898[19]) 


(0.941[19]) 



Table 3 

Ratios of anti-strange particle to strange particle in relativistic heavy ion collisions. 

Results at full RHIC oncrgy arc ol)tainod by extrapolation. 





AGS* 


SPS 


RHIC 


full RHIC 


A/A 


1.40x10-3 


0.134 


0.738 


0.930 




(2.0 ± 0.8) X 10-3 [27,9] 


0.131±0.017[24,12] 


0.77±0.07[26,19] 






1.74x10-3 


0.211 


0.823 


1.00 






0.232zb0.033[25,12] 


0.82±0.08[19] 






8.12x10-^ 


0.378 


0.926 


1.04 






0.383±0.081[24,12] 


(0.898[19]) 


(0.941[19]) 


K-/K+ 


4.54-1 


1.82-1 


0.800 


0.791 




(4.4 ±0.4)-i [28,9] 


(1.8±0.1)-i[24,12] 


0.88±0.07** 





* Both data and results are for Si+Au collision at 14. 6A GeV/c. 

* Average value of four experimental data shown in [19]. 
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